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SECTION  I 
INTRODUCTION 


This  report  presents  the  results  of  a  preliminary  evaluation  of 
the  Short-Period  (SP)  Norwegian  Seiamic  Array  (NORSAR),  using  seismic  data 

recorded  during  1970,  1971,  and  1972.  The  overall  objectives  of  the  NORSAR 
SP  evaluation  are: 

•  Determine  the  best  processing  methods  for  enhancing 
the  signal -to-noise  ratio  of  Eurasian  events 

•  Determine  the  array  detection  capability  for  Eurasian 
events 

•  Evaluate  the  performance  of  short-period  discriminants 
at  NORSAR 

•  In  conjunction  with  the  long -period  NORSAR  data, 
determine  the  detection  and  discrimination  capability 
of  NORSAR  for  Eurasian  events. 

Substantial  progress  has  been  made  toward  achieving  the  first  three  objectives; 
future  analysis  will  be  directed  toward  improving  the  preliminary  results  pre¬ 
sented  in  this  report  and  meeting  the  fourth  objective. 

Five  analysis  tasks  were  undertaken  in  order  to  meet  the  first 
three  objectives  stated  above: 

•  Noise  analysis 

•  Signal  analysis 

•  Array  processing  effectiveness 

•  Preliminary  detection  threshold  estimation 

•  Behavior  of  SP  discriminants 
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Results  of  each  task  are  presented  in  subsequent  sections  of  this  report. 

The  NORSAR  SP  array,  centered  about  100  km  due  north  of  Oslo, 
Norway,  consists  of  132  short-period  seismometers  and  has  an  aperture  of  about 
100  km.  The  sensors  are  grouped  into  22  six-element  subarrays;  each  subarray 
has  a  center  sensor  and  a  five -sensor  ring  and  is  about  7  km  in  diameter  (Fig¬ 
ure  I-l).  In  this  report  subarray  1  refers  to  subarray  OlA,  subarrays  2  through 
8  to  subarrays  OIB  through  07B,  and  subarrays  9  through  22  to  subarrays  OIC 
through  14C.  Within  a  subarray,  sensors  0  through  4  refer  to  sensors  in  the 
surrounding  ring,  starting  with  the  first  sensor  east  of  a  north  (0°)  azimuth, 
and  proceeding  clockwise  around  the  ring.  Sensor  5  refers  to  the  central  sensor 
of  a  subarray.  Thus,  sensor  0  of  subarray  1  is  sensor  OlAOl  in  the  official 
nomenclature,  and  sensor  5  of  subarray  22  is  sensor  14C00. 

The  results  presented  in  the  following  sections  are  based 
primarily  on  analysis  of  events  located  on  the  Eurasian  continent.  A  total  of  102 
events  have  been  analyzed;  Table  I-l  lists  the  event  parameters. 

Geographically,  the  events  are  concentrated  principally  in  the 
Northwestern  Pacific  (from  Kamchatka  to  Taiwan)  and  in  south  central  Asia  (north 
and  west  of  the  Himalayan  system).  There  are  five  events  from  the  Ural  Moun¬ 
tain  region,  three  from  continental  North  America,  two  from  .,\s  far  west  in 
Europe  as  Greece,  two  from  the  Arctic  Ocean  and  one  from  the  Aleutian  Islands. 
Ninety  percent  of  the  events  are  either  shallow  (<100  km)  focus  or  have  unknown 
depths.  Twelve  events  are  signals  from  known  test  areas,  including  six  from 
Eastern  Kazakh,  one  from  Novaya  Zemlya,  two  from  the  Ural  Mountains,  one 
from  the  Aleutian  Islands  and  two  from  Nevada. 

For  72  of  the  signals  an  associated  320-second  noise  sample  taken 
just  prior  to  the  signal  was  edited  and  processed,  these  samples  were  used  to 
obtain  the  noise  analysis  results  reported  in  Section  II  of  this  report. 
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FIGURE  I-l 


NORSAR  SHORT  PERIOD  ARRAY 


TABLE  1^1 

I 

EVENT  PARAMETERS 


Designation 


KUR/060/07N 

1 

TAI/062/o4N 

JAP/O63/O7N 

MED^63/08N 

JAP/065i/11N 

UZB/065/23N 

JAP/066/02N 

TAL/066/21N 

TUR/067/23N 

GRE/068704N 

CAL/0^)8/15N 

JAP/068/2iN 

I^UR/q69/05N 

SZE/070/04N 

KIR/070/24N 

ALI/07,1/08N 

MON/071/Q8N 

RYa/072/03N 

,  ' 

JAP/072/03N 

RYU/072/03N 

VAN/072/23N 

SAK/073/Q7N 

KM1/073/12N 

KM2/173/12N, 


_  .  .  ,  Source 

Origin  Time  Lat.  Long,  Depth  ^uUgtin  Comment 

06.  53.  32  50.  7N  154.  5E 

04.01:45  26.  8N  125.  9E 

07.48.39^  40.  7N  143.  5E 
08.  58.  29  33.*5N 

11.49.40  33.  lNil42.0E 

23.02.02  40.  8N 

I  I 

02.50.02  43.  ON  145.  8E 

I 

21.55.17  4.  3N  127.  IE 

22.  44.47  37.  5N 

04.  58.  38  ,  38.  7N 

1  '  ' 

15.35.16  {36.  8N  122.  2W 

'  i 

12.19.52  'i36.  6n  142.  IE 


03/,01/7i 
03/03^71, 

03/04/71 
03/04/71 
03/06/71 
03/06/71 
03/07/71 
03/07/71 

03/08/71  22.  44.47  37.  5N 

03/09/71  04.  58.38,  38.  7N 

1  '  ' 

,03/09/71  15.35.16  |36. 8N  122. 2W 

03/09/71  12.  19.52  '{36.  6n  142.  IE 

03/10/71  05.21.03  j45. 4N|  15i.2E 

03/11/71  04^43.37  28.  7N  103. 7E 

03/11/71  23.59.46  40,  IN  77.  IE 

03/12/71  08.  24.47  51.  5N,177. 4E 

03/1^71,  08.56.04  47. 7N  113.  9W 

03/13/71  02.37.40^  26. 8N  129.  8E 

03/13/71  02.  59.38  40.  IN  142. 4E 

03/13/71'  03.  32.  20  '  26.  8N  129.  7E 

03/13/71  23.  51.35  50.  6N  1  129.  9N 

03/14/71  06.41.10  46.  3N  142. 3E 

03/14/71  12.|15.  14  53.  9N  160. 5E 

'  i 

03/14/71  12.24.04  54.  ON  160. 6E  4.8 


MBN 


MBN 


I 
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TABLE  I- 1 
(Continued) 


1  Designation 

Date 

Origin  Tim 

e  Lat. 

Long. 

"'b 

Dept 

Source 
^  Bulletii 

^  Comment 

JAP/074/14N 

03/15/71 

o 

• 

• 

31.  IN 

141.  81 

:  4.8 

N 

P 

GRE/074/15N 

03/15/71 

15.  23. 17 

37.  2N 

24.  01 

:  4.8 

32 

p 

SIN/0  75/ 22N 

03/16/71 

22.  28.  33 

39.  IN 

75.  71 

:  4.5 

N 

p 

TAI/075/22N 

03/16/71 

22.  52.09 

23.  8N 

121.  6E 

4.6 

40 

P 

SWR/076/02N 

03/17/71 

02.  15.46 

46.  ON 

45.  4E 

4.  1 

N 

s 

ND 

KUR/077/03N 

03/18/71 

03.  13.  15 

49.  2N 

156.  3E 

5.  1 

N 

P 

KUR/077/08N 

03/18/71 

08.  11.  04 

49.  3N 

155.  3E 

3.9 

60 

s 

ND 

KAZ/081/04N 

03/22/71 

04.  32.  57 

49.  7N 

78.  2E 

5.  8 

0 

p 

URA/082/06N 

03/23/71 

06.  59.  56 

61.  3N 

56.  5E 

5.6 

0 

p 

KIR/082/09N 

03/23/71 

09.  52.  12 

41.  5N 

79.  3E 

5.  7 

N 

p 

KIR/082/20N 

03/23/71 

20.47.  17 

41.  5N 

79.  3E 

N 

P 

TSI/083/13N 

03/24/71 

13.  54.  17 

35.  5N 

98.  2E 

5.8 

13 

P 

KIR/083/20N 

03/24/71 

20.  54.  28 

41.  5N 

79.  5E 

5.  3 

18 

P 

AND/087/08N 

03/28/71 

08.  23.  19 

11.  8N 

95.  IE 

5.  5 

N 

P 

SIN/0  90 /20N 

03/31/71 

20.  00.  31 

39.  6n 

74.  8E 

4.  8 

38 

p 

TIB/093/07N 

04/03/71 

07.  34.  50 

32.  2N 

95.  IE 

5.  1 

N 

P 

TAD/094/01N 

04/04/71 

01.  35.  23 

38.  4N 

73.  3E 

4.  8 

N 

p 

KUR/095/00N 

04/05/71 

00.  00.04 

44.  IN 

147.  4E 

4.2 

N 

s 

KUR/095/17N 

04/05/71 

16.  49.08 

45.  5N 

151.  9E 

5.  1 

N 

p 

KUR/098/03N 

04/08/71 

03.  23.  50 

46.  9N 

151.  3E 

3.  8 

N 

s 

KUR/099/04N 

04/09/71 

03.  56.  07 

43.  4N 

147.  6E 

1.  8 

50 

p 

KAM/099/08N 

04/09/71 

08.  33.  20 

56.  3N 

162.  7E 

1.9 

45 

p 

KUR/099/15N 

04/09/71 

15.08.09 

44.  3N 

147.  OE  f 

>.  1 

L26 

p 

YUG/100/02N 

34/10/71 

02.  58.  05 

42.  5N 

m 

t.  6 

21 
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TABLE  I-l 
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Designation 

Date 

Origin  Time 

Lat. 

Long. 

■ 

Depth 

Source 

Bulletin 

Comment 

OKH/101/02N 

04/11/71 

02.  36.  53 

59.  ON 

157.  3E 

4.0 

N 

S 

ND 

KUR/101/08N 

04/11/71 

08.42.05 

46.  2N 

153.  OE 

4.  5 

N 

P 

FTIR/102/18N 

04/12/71 

18.06.  32 

43.  4N 

147.  3E 

4.  1 

35 

S 

IRA/102/19N1 

04/12/71 

19.03.  25 

28.  3N 

55.  6E 

6.  0 

44 

P 

SIB/105/ 18N 

04/15/71 

18.53.57 

70.  IN 

155.  OE 

3.  5 

N 

S 

ND 

KAZ/108/23N 

04/18/71 

23. 14. 19 

45.  8N 

76.  IE 

4.  5 

N 

S 

GRE/109/02N 

04/19/71 

02.43.  52 

39.  ON 

20.  5E 

5.  1 

16 

P 

TAD/110/03N 

04/20/71 

03.48.27 

38.  3N 

73.  5E 

4.  9 

130 

P 

i 

» 

KAZ/115/03N 

04/25/71 

03.  32.  58 

49.  8N 

78.  IE 

5.9 

0 

P 

i 

) 

( 

1 

raN/121/19N 

05/01/71 

18.58.11 

36.  3N 

70.  3E 

4.  5 

230 

P 

t 

1 

j 

TIB/123/00N 

05/03/71 

00.  33.22 

30.  8N 

84.  5E 

5.4 

16 

P 

j 

CAS/i35/04N 

05/15/71 

04.53.05 

38.  IN 

49.  IE 

4.  6 

N 

P 

URA/136/07N 

05/16/71 

06.  59.  33 

52.  ON 

56.  OE 

3.  8 

N 

S 

ND 

YUN/.141/02N 

05/21/71 

02.  58.  37 

26.  7N 

101.  8E 

4.9 

45 

P 

KUR/142/12N 

05/22/71 

12.59.  53 

49.  ON 

154.  OE 

3.  8 

N 

S 

ND 

KAZ/145/04N 

05/25/71 

04.02.  57 

49.  8N 

78.  2E 

5.  2 

0 

P 

i 

t 

TAD/147/00N 

05/27/71 

00. 30. 27 

38.  3N 

69.  OE 

4.  8 

36 

P 

TIB/155/14N 

06/04/71 

14. 10.46 

32.  2N 

95.  2E 

5.0 

N 

P 

SZM/155/09N 

06/04/71 

09. 10.02 

84.  6N 

108.  OE 

5.  1 

N 

P 

NEC/156/10N 

06/05/71 

10.21.28 

37.  3N 

113.  7E 

4.7 

N 

P 

i 

KAZ/157/04N 

06/06/71 

04.02.  57 

50.  ON 

77.  6E 

5.  5 

0 

P 

i 

NEV/167/14N 

06/16/71 

14.50.00 

37.  ON 

116.0W 

4,9 

0 

P 

KAZ/170/04N 

06/19/71 

04.03.  57 

50.  ON 

77.  7E 

5.  5 

0 

P 

KAZ/181/04N 

06/30/71 

03.  56.  57 

50.  ON 

79.  IW 

5.4 

0 

P 
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TABLE  I-l 
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Designation 

Date 

Origin  Time 

Lat. 

Long. 

Depth 

Source 

Bulletin 

Comment 

KIR/184/04N 

07/03/71 

04.26.22 

41.  3N 

79.  3E 

4.9 

17 

P 

KUR/185/15N 

07/04/71 

15.29.21 

43.  7N 

147.  5E 

4.4 

N 

P 

TAD/188/04N 

07/07/71 

03.  52.  53 

38.  6N 

73.  IE 

4.5 

62 

P 

NEV/189/14N 

07/08/71 

14.00.00 

37.  4N 

117.  IW 

5.  5 

0 

P 

KUR/191/02N 

07/10/71 

02.04.28 

51. ON 

153. OE 

3.6 

N 

S 

KUR/191/14N 

07/10/71 

14.28.56 

43.  5N 

147.  8E 

4.8 

55 

P 

URA/141/16N 

07/10/71 

16.59.59 

64.  2N 

55. 2E 

5.  3 

0 

P 

KUR/193/06N 

07/12/71 

06.08.15 

45.  ON 

149.  OE 

4.6 

N 

S 

KUR/199/12N 

07/18/71 

12.  32.38 

51. ON 

157. OE 

4.3 

N 

S 

ND 

KAM/201/05N 

07/20/71 

05.33.24 

56.  ON 

161. OE 

3.7 

N 

s 

WRS/202/15N 

07/21/71 

15.49.27 

52.  ON 

54.  OE 

3.8 

N 

s 

ND 

KUR/202/23N 

07/21/71 

23. 50. 14 

47.  ON 

147.  OE 

3.  8 

N 

s 

ND 

HOK/203/16N 

07/22/71 

16.04. 16 

41.  ON 

147.  OE 

3.6 

N 

s 

ND 

KUR/204/12N 

07/23/71 

12.45.12 

46.  ON 

150.  OE 

3.7 

N 

s 

ND 

KUR/205/04N 

07/24/71 

04. 19.41 

45.  ON 

147.  OE 

3.6 

N 

s 

ND 

SWR/205/11N 

07/24/71 

11. 11.42 

48.  ON 

28.  OE 

3.5 

N 

s 

KUR/206/00N 

07/25/71 

00.41.26 

50.  ON 

154.  OE 

4. 1 

N 

s 

ND 

KAM/206/01N 

07/25/71 

01.23. 19 

55.  ON 

163. OE 

4.0 

N 

s 

ND 

KAM/206/08N 

07/25/71 

08. 12.34 

53.  ON 

160. OE 

3.7 

N 

s 

KUR/206/08N 

07/25/71 

08.35.19 

44.  ON 

147. OE 

4.0 

N 

s 

SIN/207/01N 

07/26/71 

01.48.33 

39.  9N 

77.  2E 

6.0 

N 

p  • 

BLS/210/19N 

07/29/71 

19.40.10 

39.  6N 

30.  4E 

4.  5 

N 

p 

KUR/213/02N 

08/01/71 

02.06.06 

m 

156.  8E 

5.6 

20 

p 

KUR/214/02N 

08/02/71 

02.  57.24 

50.  8N 

158.  IE 

3.6 

N 

s 

ND 

KYU/214/13N 

08/02/71 

13.59.41 

33.  ON 

130. 9E 

3.6 

N 

s 

I 
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Designation 

Date 

Origin  Time 

Lat. 

Long. 

B 

Depth 

Source 

Bulletin 

Comment 

KAM/217/01N 

08/05/71 

01.05.57 

50.  ON 

156.  8E 

3.  7 

N 

S 

ND 

SIN/221/01N 

08/09/71 

01.03. 16 

42.  IN 

83. 4E 

4.2 

N 

P 

IRA/221/02N 

08/09/71 

02.54.37 

36.  2N 

52.  7E 

5.2 

27 

P 

SAK/248/18N 

09/05/71 

18.  55.  57 

48.  ON 

143.  OE 

4.  3 

N 

S 

TUR/252/15N 

09/09/71 

15. 10.03 

37.  3N 

30.  2E 

5.  3 

23 

P 

NZM/270/06N 

09/27/71 

05.  59.  55 

73.  4N 

55.  IE 

6.  4 

0 

P 

WRS/277/10N 

10/04/71 

10.00.02 

61. 6N 

47.  IE 

5.  1 

13 

P 

WRS/295/05N 

10/22/71 

05.00.00 

51.  6N 

54.  5E 

5.3 

6 

P 

CAN/310/22N 

11/06/71 

22.00.00 

51.  5N 

179.  IE 

6.8 

2 

P 

Abbreviations; 


S  SAAC/LASA  Bulletin 
P  PDE  Bulletin 
ND  NONDETECTION 
M^B  NORSAR  Magnitude 
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Data  quality  was  excellent;  about  one -half  of  the  time  all  132 
sensors  were  operational.  In  mosc  of  the  other  cases  one  subarray  (six  sensors) 
was  dead;  the  worst  data  loss  was  19  sensors.  With  the  exception  of  one  event, 
there  were  essentially  no  spikes  in  the  data.  We  do  not  have  individual  seis¬ 
mometer  response  curves,  however,  based  on  results  of  the  signal  and  noise 
analysis  it  appears  that  the  seismometers  are  reasonably  well  equalized  across 
the  array. 


SECTION  n 
NOISE  ANALYSIS 

A.  INTRODUCTION 

Seventy-two  320-second  (5  minutes,  20  seconds)  noise  samples 
spread  over  1971  were  analyzed. to  determine  the  NORSAR  short  period  noise 
characteristics;  all  noise  samples  preceded  events  and  were  selected  to  end  25 
seconds  before  the  predicted  P-wave  arrival  time. 

The  following  calculations  were  made  during  the  analysis: 

•  Spectral  content  and  variability  across  the  array  for  single  sensors 
and  subarray  beams 

•  RMS  noise  levels  for  selected  bandpasses  for  the  single  sensors, 
subarray  beams  and  the  array  beam 

•  Variability  with  time 

•  Multiple  coherences  both  within  a  subarray  and  between  subarrays. 

The  noise  data,  with  the  exception  of  level  as  a  function  of  time,  were  very  similar 
sample -to -sample.  Thus  the  selected  data  presented  in  this  report  is  representa¬ 
tive  of  the  entire  72 -sample  ensemble. 

B.  SPECTRAL  CONTENT 

The  average  single  sensor  spectrum,  average  subarray  beam  spec¬ 
trum,  and  array  beam  spectrum  for  the  noise  sample  preceding  the  Greece  event 
(9  March  1971)  are  shown  in  Figure  U-l.  It  should  be  noted  that  the  Greece  sam¬ 
ple  has  higher  than  average  noise  levels.  The  system  response  has  not  been  re¬ 
moved,  but  the  spectra  are  normalized  to  l(m;i)^/Hz  at  1.  0  Hz.  The  spectra  are 
very  simple,  with  a  peak  at  .  16  to  .  33  Hz  (6  to  3  seconds)  and  a  rapid  decrease  in 
spectral  density  with  increasing  frequency.  There  appears  to  be  a  weak  secondary 
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Noise  Power  Density  (dB  relative  to  1  (m/i)  /Hz  at  1.0  Hz 


1.0  2.0  4.0 


Frequency  (Hertz) 

FIGURE  11- 1 

AVERAGE  SINGLE  SENSOR  AVERAGE  SUBARRAY  BEAM  AND  ARRaV  BEAM 
NOISE  POWER  SPECTRA  FOR  GREECE  NOISE  SAMPLE 

II-2  ,  20 


0 


0 

0 

0 

D 

u 

0 

u 
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peak  (best  seen  on  the  array  beam  spectrum)  at  about  0.  8  Hz,  however  this  peak 
is  more  than  10  dB,  down  from  the  low-frequency  peak. 

i  ‘ 

,  The  single  sensor  power  spectra  variations  are  summarized  in 
Figure  lt-2,  which,  shows  a  typical  sensor  to  average  single  sensor  spectral 
ratio,  the  maximum  positive  ratio  (sensor  3  ^of  subarray  10)  and  the  maximum 
negative  ratio  (sensor  2  of  subarray  11).  Figure  lI-2  shows  that  the  132  single 
sensor  power  spectra  are  very  similar;  the  maximum  positive  or  negative  varia¬ 
tion  is  about  5  dB,  apd  for  most  sensors  it  is  less  than  2  dB  at  all  frequencies. 

'  ■  i 

,  Only  one  anomaly  was  noted  in  the  single  sensor  power  spectra; 

sensors  O.and  1  of  subarray  22  had  a  spectral  line  at  2.  8  Hz.  This  anomaly  was 
observed  in  most  of  the  remaining  noise  samples 

Variation?  in  the  noise  spectra  of  the  subarray  beams  were  ex¬ 
amined  by  comparihg  the  ratio  of  each  subarray  beam  to  the  average  subarray 
beam.  In  the  region  of  significant  poWer  density  (0  to  2.  0  Hz)  all  subarrays  were 
within4  dB  of  the  average  subarray  spectrum  (Figure  II-3).  At  higher  frequencies 

'  variation?  sometime?  we^e  larger,  apparently  because  subarray  system  noise 
levels  aie  different. 

In  summary,  the  NOR5AR  short -period  noise  spectra  are  very  sim- 
I  pie  ,aed  very  similar  (except  ii,  level)  from  sample  to  sample.  No  significant  varia- 
tions  either  in  spectral  shape  or  across  the  array  have  been  observed. 

I 

C.  RMS  NOISE  LEVELS  AND  TIME  VARIABILITY 

1  I 

Wideband  RMS  noise  levels  for  the  single  sensors,  subarray  beams 
and  array  beam  were  computed  for  all  noise  samples.  In  addition,  for  the  ref¬ 
erence*  subarray  apd  adjusted  delay  array  beams**  RMS  calculations  were  made 

*  ■  The  reference  subarray  beam  is  the  subarray  beam  selected  as  reference 

for  the  associa,ted  event;  it  i?  the  subarray  wi^h  the  largest  signal  ampli¬ 
tudes.  I 

**  !  The  adjusted-delay  arra^  beam  is  the  array  beam  formed  with  time  delay 

anomalies  appropriate  for  the  associated  event. 
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Spectral  Ratio  (dB)  Single  Sensor/ Average  Single  Ser^'or 


Frequency  (Hertz) 

FIGURE  II-2 

SINGLE  SENSOR  NOISE  SPECTRAL  RATIO  VARIATIONS 
GREECE  NOISE  SAMPLE 
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Spectral  Ratio  (dB)  Subarray  Beam/ Average  Subarray  Beam 


0  1  2  3  4  s 


Frequency  (Hz) 

FIGURE  II -3 

SUBARRAY  BEAM  NOISE  SPECTRAL  RATIO  VARIATIONS 
GREECE  NOISE  SAMPLE 
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after  application  of  two  bandpass  filters  (Figure  n.4):the  "standard”  filter,  which 
IS  a  filter  that  appears  to  be  optimum  for  improving  the  slgnal-to-noli©  ratio 
on  the  array  beam;  and  the  SDL  filter,  which  has  been  used  at  several  other 
short-period  arrays.  Thus  the  "standard"  filter  RMS  noise  levels  represent 
minimum  levels  that  can  be  achieved  (without  significantly  attenuating  the  signal). 
The  SDL  noise  levels  can  be  compared  with  RMS  values  obtained  at  other  short- 
period  arrays  such  as  LASA.  It  should  be  noted  that  the  system  response  has 
not  been  removed  prior  to  calculating  an  RMS  level,  but  that  the  levels  have  been 
normalized  to  1.  0  Hz.  using  the  42  x  lo'^  m^/  count  conversion  factor. 

Figure  II-5  is  a  histogram  of  wide-band  single  sensor  RMS  noise 
levels  for  three  noise  samples;  two  (TAD/188/04N  and  KM2/073/12N)  were  taken 
during  relatively  quiet  periods,  one  (KUR/193/06N)  was  taken  during  a  relatively 
noisy  period.  For  the  quiet  noise  samples  the  mode  is  about  6  m^,  a  few  values 
as  low  as  3  and  a  few  as  high  as  10  were  observed.  The  total  spread  is  ±  6  dB, 
which  is  consistent  with  the  spectral  variations  shown  previously.  For  the  noisy 
sample  the  mode  is  about  12  m^;  values  range  from  5  to  17  (close  to  ±  6  dB),  how¬ 
ever  the  histogram  is  less  peaked.  Of  interest  is  whether  the  low  (or  high)  values 
are  associated  with  the  same  sensors  for  all  samples  and  if  so,  if  the  low  (or  high) 
noise  values  correspond  to  low  (or  high)  signal  amplitudes.  We  will  be  studying 
this  aspect  in  the  future. 

Figures  II-6  and  II-7  show  the  wide-band,  SDL-filtered  and  standard- 
filtered  RMS  values  for  the  reference  subarray  beam  and  the  adjusted-delay  beam 
respectively  for  all  72  noise  samples.  For  the  reference  subarray,  which  has 
noise  typical  of  all  subarrays,  wide-band  RMS  noise  levels  range  from  about  0.8 
to  6  mfj,  and  typically  are  1.5  There  is  evidence  of  a  seasonal  dependence 

in  the  levels  as  summertime  levels  (days  130  to  270)  are  somewhat  lower  than 
wintertime  levels.  Analysis  of  the  long -period  NORSAR  data  (Texas  Instruments, 
1972)  showed  that  the  3  to  6  second  microseismic  peak,  which  dominates  these 
wide-band  RMS  measurements,  increases  significant.,  uuring  the  winter;  thus  this 
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FIGURE  U-4 

FREQUENCY  RESPONSE  OF  STANDARD  FILTER 
AND  SDL  FILTER 
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KUR/193/06N 


Si 


I 


sur  irtt^iiaiBi^i 


wlif  hind 
SDh  FUter 
SLindird  Filter 
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Day  of  Year  (r971)  , 

FIGURE  II-6(a)  i 

REFERENCE  SUBARRAY  NOISE  LEVELS 


170  180  190  ,  206 


RMS  (m/x) 


42.0, 


4.2. 


200 
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Day  of  Year  (1971) 

FIGURE  n-6(b) 

REFERENCE  SUBARRAY  RMS  NOISE  LEVELS 


RMS  (m/i) 


1.2ft 


.04 
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Day  of  Year  (1971) 

FIGURE  n-7(a) 

ADJUSTED -DELAY  ARRAY  BEAM  RMS  NOISE  LEVELS 
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seasonal  dependence  would  be  expected.  The  SDL  filter,  with  its  low  corner 
frequency,  does  not  reduce  the  reference  subarray  beam  RMS  noise  levels  very 
much.  RMS  levels  reuige  from  0.  6  to  3  mfj,  and  typically  are  1  mfji.  The 
"standard"  filter,  which  has  a  very  severe  low-frequency  roll-off,  attenuates 
the  subarray  beam  RMS  noise  levels  significantly,  as  would  be  expected.  Values 
range  from  0.  2  to  0.  8  mfj,  and  typically  are  0.  5  mfi.  Thus  this  filter  reduces 
RMS  noise  levels  by  10  dB,  and  more  on  noisy  days. 

For  the  adjusted-delay  array  beam  wide -band  RMS  noise  levels 
range  from  0.25  to  1.2  mjj,  and  typically  are  0.5  mj[|.  Again,  a  seasonal  depen¬ 
dence  is  evident.  Note  however,  that  there  are  occasional  sumer  days  (e.  g.  ^ 
day  193)  when  the  noise  levels  are  high.  RMS  levels  through  the  SDL  filter  range 
from  0.  16  to  0.  6  mfji,  and  typically  are  0.  3  mfj.  RMS  levels  through  the  "standard" 
filter  range  from  0.  04  to  0.  17  mfji,  and  typically  are  0.  12  mfji.  The  "standard" 
filter  is  probably  close  to  the  best  filter  for  detection  purposes  (see  Section  V)  and 
is  similar  to  the  filter  currently  being  used  on  the  NORSAR  on-line  system.  RMS 
levels  through  this  filter  are  about  a  factor  of  2  higher  than  those  observed  at  LASA 
(Dean,  1971)  which  implies  the  NORSAR  teleseismic  detection  threshold  will  be 
about  0.  3  m^^  units  higher  than  the  3.  9  (90%  incremental)  LASA  threshold. 

Finally,  it  is  interesting  to  note  that  the  seasonal  dependence  is 
less  evident  on  the  filtered  beams;  this  implies  that  the  low-frequency  microseis- 
mic  noise  peak  shifts  to  longer  periods  when  the  levels  increase. 

D.  MULTIPLE  COHERENCE 

Noise  multiple  coherences  within  a  subarray  (five  sensors  on  ring 
predicting  the  central  sensor  of  the  subarray)  were  calculated  for  the  22  sub¬ 
arrays  for  several  noise  samples.  Atypical  multiple  coherence  plot  (subarray  10) 
is  given  in  Figure  U-8  for  the  Greece  noise  sample.  Except  for  the  peak  at  3  to  6 
seconds,  coherence  is  low  at  aU  frequencies  where  there  is  significant  energy.  The 
3  to  6  second  energy  is  well  below  the  signal  band,  and  so  can  be  removed  by  band¬ 
pass  filtering.  Above  2  Hz,  there  is  moderate  coherence,  but  spectral  levels  are 
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MULTIPLE  COHERENCE  FOR  SUBARRAY  10  NOISE  FIELD 
GREECE  NOISE  SAMPLE 
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at  least  25  dB  below  the  peak.  The  moderate  coherence  at  high  frequencies  may 
be  due  to  the  system  electronics.  In  any  event  it  is  clear  that  multichannel  pro¬ 
cessing  would  not  be  necessary  at  the  subarray  level. 

Multiple  coherences  calculated  between  subarrays  were  very  low, 
as  expected,  because  of  the  large  distances  involved.  Figure  II-9  shows  the 
multiple  coherence  between  the  reference  sensor  of  subarray  1  and  the  reference 
sensors  of  subarrays  2,  3,  4,  6,  and  8  for  the  Greece  noise  sample;  very  low 
values  were  obtained  over  the  entire  band. 
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FIGURE  II- 9 

MULTIPLE  COHERENCE  BETWEEN  THE  REFERENCE 
SENSOR  OF  SUBARRAY  1  AND  THE  REFERENCE 
SENSORS  OF  SUBARRAYS  2,  3,  4,  6  AND  8 
GREECE  NOISE  SAMPLE 
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SECTION  m  I 

SIGNAL  ANALYSIS 

I 

A.  INTRODUCTION 

Several  large  signals  recorded  by  the  entire  array  were  analyzed 
in  order  to  study  signal  characteristics.;  The  following  lanalyses  were  iperfbrmed:  • 

•  Qualitative  review  of  signal  similarity  among  sensors  at  a  subarray 

and  eimong  subarrays  , 

•  Subarray  beam  amplitude  variations  '  ! 

i 

•  Spectral  content  ' 

•  Time  delay  anomalies  (deviation  from  plane  wave  propagation) 

•  Comparison  of  NORSAR  magnitudes  (m^^)  with  PDE  (or  LAS  A) 
magnitudes 

Results  of  these  analyses  are  discussed  in  subsequent  subsections  ‘  i 

B,  SIGNAL  SIMILARITY 

I  I 

'  .  '  i 

Qualitative  estimates  of  the  fiingle  sensor  signal  similarity  within  a 
subarray  were  made  by  comparing  the  single  sensor  waveforms  of  the  Greece  and 
Vancouver  Island  events  for  all  22  subarrays.  In  general,  waveforms  were  quite 
similar  within  a  subarray  (In  Section  IV  data  are  presented  which  shows  that  the 
signal  loss  in  subarray  beamforming  typically  was ,  1  dB;  this  small  loss  is  consis¬ 
tent  with  the  above  qualitative  observations).  Some  nearby  (A  <30°)  Eurasian 
events  contained  significant  amounts  of  high-frequency  energy  (f  >2.  0  Hz)  which 
was  not  preserved  by  the  subarray  beam,  probably  because  the  time  delay  anomaly 
corrections  were  not  precise  enough  for  these  events. 

Signal  similarity  among  subarrays  was  quite  variable  from  event 
to  event.  For  predominantly  low-frequency  events  (Figure  III-l)  similarity  was' 
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generaUy  good  over  the  entire  wavetrain.  For  events  with  substantial  high- 
f^quency  energy  (Figure  IH-2).  however,  subarray  wavetrains  often  were  quite 
different  after  the  first  few  cycles  (usually  the  low-frequency  part  of  these  sig¬ 
nals  was  not  significantly  attenuated  in  array  beamforming,  implying  that  signal 
similarity  at  low  frequencies  was  reasonably  good).  A  detailed  discussion  of 
array  beamforming  loss  is  given  in  Section  IV. 

C.  SUBARRAY  BEAM  AMPLITUDE  VARIATIONS 


Signal  amplitude  variations  among  subarray  beams  are  striking 
(note  that  the  seismometers  appear  to  be  well-equalized).  For  the  GRE/068/04N 
event,  the  peak-to-peak/2  values  are  shown  in  Figure  111-3.  These  values  vary 
from  84  m^at  subarray  12  to  24  m„  at  subarray  6  (10.8  dB  variation).  For  this 
event,  the  southwest  quadrant  of  the  NORSAR  array  (subarrays  6.  7, 15,  16,  17, 18, 
19.  20.  and  21)  exhibits  peaks  below  42  m^  .  whUe  the  northeast  section  (sub¬ 


arrays  12.  13,  and  14)  has  peaks  above  63  mp.  For  the  URA/082/06N  event,  the 
corresponding  zero-to-peak  values  are  shown  in  Figure  111-4.  These  values  vary 
from  530  at  subarray  11  to  228  m,i  at  subarray  7  (3.  2  dB  variation).  Sub¬ 
array  12  also  had  the  highest  signal  amplitude  for  this  event. 


Table  m-l  shows  the  ratio  (in  dB)  of  the  strongest  to  weakest  sub¬ 
array  RMS  signal  level  for  events  having  m^  >  5.  0.  The  signal  RMS  calculations 
were  made  over  a  6. 4  second  gate  which  began  just  prior  to  the  P-wave  arrival. 
The  events  are  grouped  by  region,  and  the  subarrays  having  the  strongest  and 
weakest  signal  respeccively  also  are  listed. 


The  average  spread  in  signal  amplitudes  was  II.  8  dB  (a  factor  of  4), 
however  the  spread  shows  strong  regional  dependence.  The  largest  spreads  were 
for  the  Kazakh  events  (about  10:1),  the  lowest  for  the  Ural  events  (about  2:1);  events 
from  both  these  areas  were  presumed  explosions. 

Note  also  that  within  a  region  the  spread  varies  considerably;  this 
is  true  even  for  the  Kazakh  events  which  are  all  located  very  close  together. 
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SUBARRAY  11  REFERENCE 


SUBARRAY  13  BEAM 
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SUBARRAY  20  BEAM 


SUBARRAY  20  REFERENCE 

FIGURE  UI-2  I - 5  sec - 1 

SET  "ECTED  SUBARP.AY  BEAM  AND  SINGLE  SENSOR  TRACES 

OF  EVENT  KAZ/115/03N 
SCALE  FACTOR  =  0.  00005  INCHES/COUNT 
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FiGURElilI-3 

SUB  ARRAY  SIGNAL  AMPLITUDES  FOR  GRE/068/04N 
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I  '  TABLE  Illi-l  ' 

I  i  '  '  ' 

MAXIMUM  VARIATION  OF  SUBARRAY  BEAM  SIGNAL 
RMS  LEVELS  ACROSS  NORSAR  ' 

!  '  (PAGE  1  OF  2)  ' 


Event 

1  1 

A® 

1 

1 

Azimuth° 

RMS 

Variation 
Max/ 'Min 
(dB) 

_ ^ _ _ _ 

Maximum 
!  Subarray 

1 

Minimum 

Suibarray 

RUR/099/15N 

1 

5.  1 

69.  0 

32.  1 

» 

18.2 

1 

3 

15  1 

KUR/O77/03N 

'  5.  1 

66.  5 

23.8 

9.8 

9 

15 

KUR/0^5/17N 

5.,1 

69.  1  ^ 

■  28. 1 

6.9 

1  3  I  ' 

21 

KUR/213/02N 

5.6 

65.5 

23.1 

12.4 

3 

19 

JAP/063/07N 

5.2 

71.4 

36.  o' 

11.5 

9 

1^> 

JAP/  072/ 03, N 

5.  1 

71.  6 

37.1 

10.6 

9  : 

‘  >  i 

6 

KMI/073/12N 

! 

5.  3 

62.8 

19.6 

14.5' 

3 

15 

KIR/082/09N 

5.7 

44.  6 

'  8.3. 2 

'  10.4 

12 

'  » 

7 

KIR/082/20N 

6.0 

44.  6 

83.2 

9.3 

12 

'  7 

KIR/ 08  3/ 2  ON 

5.  3! 

44.  7 

83.0: 

15.  1 

11 

7  1 

SIN/20i;01N 

1 

■  6. 0 

44.8 

'86.4 

9.7 

11 

7 

KAZ/181/04N 

5.4 

38.  3 

74.5 

1  1 

16. 1 

14 

20 

KAZ/081/04N 

5.8 

38.  1 

75.5 

20.7 

14 

20,  ' 

KAZ/170/04N 

5.5 

37.  6 

75.5 

'  25.3 

14 

I 

20 

KAZ/115/03N 

5.^ 

38.0 

75.4 

23.  3 

'll  ' 

20 

KAZ/145  /  04N 

5.2 

38.  0 

75.4 

20.6 

14 

20! 

KAZ/157/04N- 

5.5 

1  ' 

'  37.7  . 

,  75.4 

l9.7 

14  ' 

20 

SZE/070/04'n 

5.  1 

66.  6 

72>  7  1 

4.8 

10 

12  1 

TIB/093/07N 

5.  1 

59.6 

77.5 

5.  9 

9 

'20 

TIB/123/00N 

5$.  6 

87.3  ' 

9.7 

7 

TIB/155/14N 

5.  0 

59.6 

.  77.5 

9 

7 

TSI/083/13N  1 

5.8 

'58.  3 

72.  9 

:  8.1 

10 

1 

7 

IRA/102/19N1 

6.  0' 

'  44.1, 

110.9 

8.8 

18 

9  ' 

IRA/i221/p2N 

5.2 

36.  0  , 

113.7 

i 

11.9, 

4 

Zz 

URA/082/06N 

5.  6 

21.6 

68:6 

5:2 

12  1 

7 

URA/191/16N 

i 

5.3 

20.  3 

61.2 

6.5 

9 

lOi  i 
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i! 

11 
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il 

u 

1! 

,0 
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TABLE  m-1 


MAXIMUM  VARIATION  OF  SUBARRAY  BEAM  SIGNAL 
RMS  LEVELS  ACROSS  NORSAR 
(PAGE  2  OF  2) 


Event 

"'b 

A° 

.Azimuth^ 

RMS 

Variation 
Max/ Min 
(dB) 

Maximum 

Subarray 

Minimum 

Subarray 

WRS/295/05N 

5.3 

25.4 

91.4 

9.4 

22 

7 

WRS/277/10N 

5.  1 

17.  3 

71.5 

11.3 

5 

7 

GRE/109/02N 

5.  1 

22.6 

160.2 

11. 1 

3 

20 

TUR/252/15N 

5.3 

26.5 

143.8 

9.2 

8 

13 

VAN/  072/ 2  3N 

5.7 

64.2 

333.5 

10.2 

1 

13 

Average  Variation  (All  Events)  =  11.8 
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The  subarrays  having  the  largest  and  smallest  signal  amplitudes 
tend  to  repeat  for  a  given  region,  especiaUy  if  general  rather  than  detailed  am¬ 
plitude  values  are  considered.  That  is.  for  the  Kuriles  -  Kamchatka  -  Japan 
region,  subarrays  3  and  9  had  high  amplitudes  for  al.’  the  events  and  subarrays 
6,  15,  16,  19.  and  21  had  low  amplitudes.  Thus  diversity-stack  weights  (for 
array  beamlorming)  calculated  from  large  events  should  be  satisfactory  for  app 
lication  to  small  events,  although  they  probably  won't  be  exactly  correct. 

FinaUy,  it  is  interesting  to  note  that  the  largest-amplitude  sub¬ 
arrays  for  some  events  are  the  smallest-amplitude  for  others.  For  example, 
subarray  9  has  large  signal  amplitudes  for  the  Kuriles  -  Kamchatka  -  Japan  re- 
gion,  and  small  signal  amplitudes  for  Iran. 

D.  SIGNAL  SPECTRAL  CONTENT 


Signal  spectral  content  was  studied  by  computing  power  density 
spectra  over  a  6.4  second  gate  beginning  just  before  the  P-wave  arrival  for  sev¬ 
eral  large  signal-to-noise  ratio  events.  Figures  III-5  to  III-8  show  the  reference 
subarray  and  adjusted  delay  array  beam  for  two  earthquakes  and  two  presumed 

explosions.  The  system  response  has  not  been  removed,  but  the  spectra  have  to 
be  normalized  to  dB  re  /Hz  at  1.  0  Hz. 


Figures  m-5  and  m-6  are  spectra  for  two  presumed  explosions, 
URA/082/06N  (A  =  22  )  and  KAZ/115/03N  (A  =  36°).  The  spectra  are  essentiaUy 
flat  out  to  about  2.  5  Hr.  then  drop  off  rapidly  at  higher  frequencies.  The  lower 
spectral  values  on  the  adjusted-delay  beam  are  not  because  of  signal  loss,  since 
the  reference  subarray  beam  is  Uie  subarray  with  the  largest  signal  amplitude. 

The  reference  subarray  signal  RMS  levels  are  3  to  13  dB,  and  typicaUy  6  dB  above 
Uie  average  subarray  level.  Thus  the  average  subarray  beam  spectra  for  these 
bvo  events  would  be  about  6  dB  lower  than  the  reference  subarray  spectra  and  hence 
It  can  be  seen  that  there  was  little  signal  loss  out  to  2.5  Hz  for  both  URA/082/06N 
and  KAZ/115/03N.  It  also  should  be  pointed  out  that  some  Eurasian  earthquakes 
especially  those  with  A  <  30°,  have  significant  high-frequency  energy. 
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Figures  III-7  and  111-8  are  spectra  for  two  earthquakes,  KIR/082/20N 
(A  =  45°)  and  TSI/083/13N  ( A  =  58°).  Both  events  have  spectral  peaks  at  about  0,8 
Hz;  KIR/082/20N  has  a  secondary  peak  at  1.4  Hz  whereas  TSI/083/13N  drops  off 
uniformly  at  higher  frequencies.  Again,  array  beamforming  signal  loss  appears  j 
to  be  small  for  both  events. 

There  appears  to  be  a  significant  difference  in  signal  energy  con-  | 

I 

tent  for  events  from  the  western  hemisphere  as  compared  with  those  from  Eurasia. 

( 

For  our  limited  ensemble  of  western  hemisphere  events,  spectral  content  generally  I 

! 

is  of  a  lower  frequency  for  both  earthquakes  and  presumed  explosions.  Figure  III-9  | 
shows  the  spectrum  for  an  event  from  the  Nevada  test  site  (A  =  73°);  this  presumed 
explosion  spectrum  drops  off  at  about  1.4  Hz.  j 

We  do  not  have  enough  data  as  yet  to  make  comments  about  spectral  | 
content  as  a  function  of  region  in  Eurasia,  except  to  say  that,  not  unexpectedly, 
Eurasian  signals  (both  earthquakes  and  presumed  explosions)  from  A  <30°  gener¬ 
ally  have  higher  frequency  content.  Regional  characteristics  of  signal  spectra 
will  be  studied  in  the  future. 

E.  TIME  DELAY  ANOMALIES 

It  was  established  early  in  the  analysis  that  standard  plane -wave 
delays  generally  were  adequate  for  subarray  beamforming,  but  that  time  delay 
anomalies  (deviations  from  plane  wave  propagation  along  the  great  circle  azimuth) 
had  to  be  accounted  for  in  array  beamforming. 

Inter -subarray  time  delay  anomalies  were  calculated  by  computing 
the  crosscorrelation  functions  between  the  reference  subarray  beam  and  the  re¬ 
maining  subarray  beams  for  all  large  SNR  events.  A  short  signal  gate  (generally 
6.4  seconds)  was  used  in  the  computations.  In  cases  where  inconsistent  or  un¬ 
certain  delays  were  obtained,  results  were  hand  checked  by  measuring  delay  an¬ 
omalies  from  the  time  traces.  After  computation,  the  delay  anomalies  were  ad¬ 
justed  so  that  they  were  all  relative  to  subarray  1  at  the  center  of  NORSAR. 
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Pow»r  Density  (dB  relative  to  l(mju)^/Hz  at  1.  0  Hz) 
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FIGURE  m-7 

REFERENCE  SUBARRAY  AND  ARRAY  BEAM  SIGNAL  SPECTRA 

OF  EVENT  KIR/082/20N 
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Power  Density  (dB  relative  to  l(in/|)  /Hz  at  1,0  Hz) 


FIGURE  in-8 

REFERENCE  SUBARRAY  AND  ARRAY  BEAM  SIGNAL  SPECTRA 

OF  EVENT  TSI/ 083/1 3N 
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Power  Density  (dB  relative  to  Hm^)  /Hz  at  1.  0  Hz) 


FIGURE  m-9 


REFERENCE  SUBARRAY  AND  ARRAY  BEAM  SIGNAL  SPECTRA 
!  OF  EVENT  NEV/ 189/ 14N 
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For  each  subarray  the  delay  anomalies  obtained  for  all  events 
were  plotted  on  a  map  of  the  Eurasian  continent  at  the  geographic  locations  (epi¬ 
centers)  of  the  events.  The  anomalies  then  were  reviewed  to  determine  the  ab¬ 
solute  value  of  the  anomalies  and  their  consistency  for  events  in  the  same  region. 
The  following  was  observed: 

•  Anomalies  typically  were  (+)  a  few  tenths  but  occasionally 
reached  a  full  second.  Anomalies  generally  were  larger 
for  subarrays  on  the  outer  ring. 

•  Consistent  delays  could  be  obtained  for  Kamchatka/ 

Kuriles,  Eastern  Kazakh,  South  Central  Russia,  and  China. 

Table  111-2  shows  the  anomalies  for  four  large  South 
Central  Russia  events;  the  individual  event  values  differ 
from  the  average  value  by  more  than  0.  1  seconds  only 

for  subarray  7.  Because  the  absolute  value  of  delay 
anomalies  overall  other  subarrays  is  0.7  seconds  or  less, 
it  appears  that  the  subarray  7  delays  determined  for  the 
Kirgiz  events  are  erroneous.  Delay  anomalies  obtained 
from  other  smaller  events  more  widely  scattered  through¬ 
out  the  area  agreed  closely  with  those  shown  in  Table  111-2. 

It  is  interesting  to  note  that  all  anomalies  were  zero  or 
negative;  this  implies  that  the  signal  was  early  (relative  to 
all  other  subarrays)  at  subarray  1. 

•  Consistent  delays  could  not  be  obtained  from  close-in  events 
from  the  Mediterranean  and  Western  Russia.  Table  111-3 
shows  the  delays  for  four  moderate  Mediterranean  events 
which  are  located  in  a  circle  with  a  radius  of  5°.  Delay 
anomalies  vary  by  as  much  as  1.  0  second  at  some  sub¬ 
arrays  event  to  event.  Anomalies  are  not  consistent  even 
for  GRE/109/02N  and  GRE/068/04N,  which  have  essen¬ 
tially  the  same  epicenter.  The  difficulty  is  obtaining  con¬ 
sistent  delays  for  the  close-in  events  results  because  the 
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events  are  rich  in  high  frequencies.  These  higher  fre¬ 
quency  arrivals  may  have  different  paths  (perhaps  be¬ 
cause  of  local  crustal  effects)  or  be  incoherent  among 
subarrays;  either  case  would  make  obtaining  consistent 
delays  very  difficult.  More  work  is  netded  to  determine 
appropriate  delays  for  these  close-in  events. 

F .  COMPARISON  OF  NORSAR  m  VALUES  WITH  PDE  AND  LASA  m  's 

D  b 

NORSAR  m^  values  were  measured  for  the  85  detected  events 

using  the  formula: 

=  log  -^  +  B 

where:  A  is  the  maximum  peak-to-peak  signal  amplitude  in  m^  on 

the  adjusted  delay  array  beam 
T  is  the  period  of  the  cycle  with  the  meiximum  amplitude 
B  is  the  distance  factor 

Values  for  B  are  shown  in  Table  II1-4;  these  also  are  used  for  calculation  of 
LASA  m,  's, 

D 

Figure  II-IO  is  a  plot  of  NORSAR  ni^'s  versus  either  PDE 
(circles)  or  LASAm^'s  (crosses).  Note  that  on  the  average  LASA  and  PDE  ni^'s 
are  the  same  (Dean,  1971)  so  that  the  data  points  can  be  considered  together. 

Values  for  presumed  explosions  are  circled. 

In  general  the  NORSAR  m, 's  are  somewhat  smaller  than  the  LASA 

D 

and  PDE  nij^'s;  the  difference  is  somewhat  more  pronounced  for  m^  <  5.  3.  We  do 
not  have  an  explanation  for  the  apparently  larger  difference  at  lower  magnitudes; 
it  is  neither  distance-  nor  region-dependent.  Further,  our  data  do  not  indicate  that 
beamforming  loss  is  lower  for  large  events.  It  may  be  simply  statistical  scatter, 
or  it  may  be  related  to  the  fact  that  several  of  the  larger  events  were  presumed 
explosions.  The  most  likely  explanation  is  that  the  PDE  magnitudes  for  small 
events  are  based  on  the  quieter  stations  in  the  PDE  net  (e.  g. ,  TFO),  which  tend  to 
give  lower  than  average  magnitudes.  We  plan  to  investigate  this  difference  in  more 
detail  in  the  future. 
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TABLE  III-4  . 

DISTANCE  FACTOR  (B)  FOR  COMPUTING  NORSAR  VALUES 

b 


Distance 

Degrees 

B 

Distance 

Degrees 

B 

1 

Distance 

1 

Degrees 

■'  B, 

0 

1 

0.  50 

35 

3.  34 

69 

3.48 

2 

1.  50 

36 

3.  34  : 

3. 49 

3 

2.  50 

37 

3.  34  , 

71  1 

3.50 

4 

2.  55 

38 

3.  33 

72 

1  3. 5o 

5 

2.  76 

39 

3.33 

73 

3.  51 

6 

2.  90 

40 

3.  32 

'  74 

3'.  51  , 

7 

3.02 

41 

3.  3i 

75 

3.52 

8 

3.  10 

42 

'  3.  32 

76 

3.53 

9 

3.  15 

43 

3.  33 

77 

3153 

10 

3.  20 

44 

3.  33  ’ 

'  78 

3.  54 

11 

3.  23 

45  . 

3.  34 

79 

.3.54  , 

12 

3.  25 

46 

3.  34 

80 

3.  55 

13 

3.  26 

47 

3.35  ' 

8^ 

3.  56 

14 

3.  26 

48 

3.  36 

82 

'  3.57 

15 

3.  25 

49 

3.  36 

83 

3.58 

16 

3.  21 

50 

'  3. 37 

84' 

3.'59  ' 

17 

3.  10 

51 

3.37  .. 

85 

,  .  3.61 

18 

2.  98 

52 

3.  38 

86 

3.64 

19 

2.  79 

53 

3.39 

87 

^.66 

20 

2.  77 

54 

3.'39 

88 

3.68 

21 

2.  80 

55 

3.40 

89  ' 

3.72 

22 

2.  85 

56 

3.40  ' 

90 

3.76 

23 

2.  94 

57 

3.  41 

91 

,3.  80 

24 

3.  04 

58 

3.  42 

92 

3.85 

25 

3.  15 

59 

3.  42 

93 

3.90 

26 

3.  24 

60 

3.  43 

94 

3.96 

27 

3.  34 

61 

3.44 

95  ^ 

4.02 

28 

3.  42 

62 

3.44  ■ 

96 

4.  .11 

29 

3.44 

63 

’  3. 45 

97 

4. 19  , 

30 

3.  42 

64 

3.45 

98 

4.28 

31 

3.  38 

65 

3  .  46 

99 

'  ■  4.  40 

32 

3.  36 

66 

3.46 

100 

4.56 

33 

3.  36 

67 

3.  47 

} 

1 

34 

3.  35 

68 

3.  48 

' 

1  .  '■ 

III-20 

52 


iiiliyiiiiiyigisg 


lassiii 


®  Pniumsd  Ejiplopion 


Figure  III- 11  is  a  histogram  of  the  magnitude  difference  (NORSAR 
minus  PDE  or  LASA).  The  shaded  area  of  the  histogram  corresponds  to  LASA 
m^'s.  There  is  a  0.  2  magnitude  unit  negative  bias;  this  could  be  accounted  for 
by  signal  loss  in  beamforming,  which  typically  is  3  to  4  dB  (0.  15  to  0.  2  magni¬ 
tude  units).  Note  that  the  LASA  magnitudes  are  corrected  for  beamforming  loss 
but  our  NORSAR  magnitudes  are  not.  The  distribution  about  this  bias  appears  to 
be  approximately  normal  and  reflects  normal  variance  in  signal  amplitude  due  to 
radiation  patterns  and/or  propagation  effects.  It  appears  that,  if  signal  loss  in 
beamforming  is  accounted  for,  NORSAR  m^'s  would  be  equivalent  to  PDE  or  LASA 
Still  to  be  explained,  however,  is  the  apparent  magnitude -dependence  of 

the  NORSAR/PDE  (or  LASA)  m^  difference. 

b 


FIGURE  m-11 

HISTOGRAM  OF  DIFFERENCES  BETWEEN  NORSAR  AND 

PDE  OR  LASA  FOR  85  EVENTS 
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SECTION  IV 

ARRAY  BEAMFORMING  RESULTS 

A.  INTRODUCTION 

This  section  covers  signal -to-noise  ratio  improvement  achieved 
by  NORSAR  short-period  array  beamforming.  Subsection  B  covers  noise  re¬ 
duction  both  at  the  subarray  level  and  at  the  array  level,  while  subsection  C 
analyzes  signal  degradation  at  both  these  levels.  Subsection  D  presents  direct 
signal-to-noise  ratio  measurements  for  subarray  beams  and  the  various  types  of 
array  beams  and  combines  these  with  the  results  of  the  two  previous  subsections 
in  an  attempt  to  specify  where  noise  redaction  and  signal  degradation  occur.  In 
addition,  the  signal-to-noise  ratio  on  the  subarray  with  the  strongest  signal  arrival 
is  compared  with  that  of  the  adjusted-delay  beam.  Finally,  the  amount  of  signal- 

to-noise  ratio  improvement  obtained  by  the  diversity- stack  beamforming  technique 
is  examined. 

B.  NOISE  REDUCTION 

Noise  reduction  in  forming  subarray  beams  and  the  array  beam  at 

NORSAR  was  only  slightly  below  the  Vn"  reduction  expected  from  sun- ming  N 
uncorrelated,  equal- strength  traces. 

Figure  IV-1  shows  noise  reduction  from  single  sensor  to  subarray 
as  a  function  of  frequency  for  typical  five-minute  noise  sample  preceding  an  event 
from  northeastern  China.  At  the  top  are  the  average  single -sensor  power  spectrum 
and  the  average  subarray  beam  power  spectrum.  The  reduction  in  noise  level  from 
single  sensor  to  subarray  beam  is  shown  at  the  bottom.  At  frequencies  above  0.  5Hz, 
it  is  remarkably  close  to  the  7.  8  dB  expected  for  six  equalized  traces. 

Figure  IV-2  gives  the  corresponding  noise  reduction  from  subarray 
to  array  beam.  At  the  top,  the  average  subarray  beam  power  spectrum  is  plotted 
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FIGURE  IV- 1 


SINGLE  SENSOR  TO  SUBARRAY  NOISE  REDUCTION 
FOR  NOISE  SAMPLE  BEFORE  EVENT  NEC/ 156/ ION 
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together  with  the  plane-wave  array  beam  spectrum  for  the  same  noise  sample. 

Noise  reduction  is  shown  at  the  bottom.  At  all  frequencies  except  0  Hz.  it  closely 

follows  the  13.  0  dB  level  expected  for  the  twenty  subarrays  operational  during  this 
sample. 

Figure  IV-3  is  a  histogram  of  the  drop  in  noise  level  from  reference 
subarray  to  adjusted-delay  beam  over  the  band  0.  55  to  1.  5  Hz  for  45  events  in 
which  all  22  subarrays  were  operational.  The  purpose  of  this  figure  is  to  indicate 
the  amount  of  noise  reduction  which  can  be  expected  from  subarray  to  array  beam. 
The  most  common  reduction  was  about  13  dB.  slightly  below  the  13.4  dB  level  ex¬ 
pected  for  22  perfectly  equalized  subarrays.  Since  one  subarray  was  used  instead 
of  all  22  subarrays,  part  of  the  spread  in  values  is  due  to  fluctuations  of  the  refer¬ 
ence  subarray  relative  to  the  average  subarray  power  in  the  band  0.55  to  1. 5  Hz. 

Noise  reduction  of  7.8  dB  at  the  subarray  level  together  with  13  dB 
at  the  array  level  yields  a  total  noise  reduction  of  slightly  below  21  dB. 

In  addition  to  demonstrating  that  noise  reduction  is  achieved 

at  NORSAR.  these  results  indicate  that  noise  levels  from  sensor  to  sensor  do  not 
vary  significantly. 

C.  SIGNAL  DEGRADATION 

To  determine  single  sensor  to  subarray  signal  degradation,  hand 
measurements  were  made  of  peak-to-peak  displacement  in  the  reference  sensor 
and  subarray  beam  for  each  subarray.  Corresponding  cycles  in  the  traces  were 
used  in  these  measurements.  Amplitude  degradation  ratios  were  averaged  across 
all  subarrays.  To  find  subarray  to  array  signal  degradation,  peak-to-peak  varia- 
t.on  was  also  measured  on  the  adjusted-delay  beam  and  compared  with  the  average 
peak-to-peak  displacement  across  the  subarray  beams.  Table  lV-1  gives  the  sub¬ 
array  level  and  array  level  signal  degradation  obtained  in  this  way  for  seven  large 
events.  The  WRS/295/05N  adjusted  delays  are  suspect.  The  subarray  beam  to 
adjusted-delay  beam  degradation  for  this  event  is  probably  the  worst  encountered 
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FIGURE  IV -3 

NOISE  REDUCTION  FROM  REFERENCE  SUBARRAY 
TO  ADJUSTED-DELAY  BEAM  (0.  55  TO  1.  5  Hz) 
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TABLE  IV-1 

I  I 

SUBARRAY-LEYEL  AND  ARRAY-LEVEL  SIGNAL  DEGRADATION 
(FROM  HAND  MEASOrSMENTS) 


, 

Event 

'Single  Sensor  to 
Subarray  Beam 
Degradation 
(dB) 

Subarray  Beam  to  • 
Adjusted-Delay 
Beam,  Degradation 
(dB) 

Single  Sensor  to 
Adjusted-Delay 
Beam  Degradation 
(dB) 

KAZ/115/03N 

,  1.4 

j 

0.  3  , 

1.7 

TIB/i23/00N 

l.'l 

1  !■ 

’  3.  1 

4.2 

TSI/083/13N 

0. 4  ’• 

.  0.3 

0.7 

VAN/ 072/ 031^ 

6.8 

0.4 

1.2 

IRA/221/02N 

1  , 

1,9 

0.3 

2.2 

WR^/295/05N 

,  1.8  ,  ’ 

11.  0 

12.8 

KAZ/081/04N 

i 

I  ' 

1.4 

1.0 

1 

2.4 

I 


I 

I 
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out  ol  all  85  signal,  detected.  Typical  signal  degradation  in  forming  subarray 
beams  appears  to  be  sUghtly  higher  than  1  dB.  Degradation  in  peak-to-peak  dis- 
placement  from  subarray  to  adjusted-delay  beam  appears  to  be  a  surprisingly 
low  1  to  3  dB  (ignoring  WRS/295/05N). 

Signal  degradation  from  subarray  beam  to  adjusted -delay  beam 
also  was  computed  using  RMS  measurements  over  a  6.4  second  gate  starting  a 
few  deciseconds  before  the  signal  arrival.  RMS  values  were  averaged  over  all 
subarray  beams  and  compared  with' the  adjusted -delay  beam  RMS.  Table  IV-Z 
shows  subarray  beam  to  adjusted-delay  beam  degradation  for  six  events  whose 
average  subarray  beam  signal -to -noise  ratio  was  greater  than  10  dB.  Results 
are  for  unfiltered  traces  and  traces  passed  through  the  "standard"  detection  filter 
(see  Figure  II-4).  Broadband  degradation  is  not  given  for  the  event  BLS/210/19N 
because  the  signal -to -noise  ratio  for  the  average  subarray  beam  was  less  than 
10  dB.  On  a  broadband  basis,  signal  degradation  was  1  to  3  dB,  except  for  the 
two  events  from  western  Russia  (which  probably  have  inaccurate  adjusted-delays). 

These  values  agree  reasonably  well  with  the  amplitude  measurements  made  directly 
from  the  time  traces. 

For  beams  passed  through  the  "standard"  detection  filter  signal 
degradation  from  subarray  to  adjusted-delay  beam  is  3  or  4  dB.  Slightly  higher 
degradation  is  observed  because  the  low  frequency  energy,  which  is  generally  more 
coherent  across  the  array,  is  attenuated  by  the  detection  filter  and  hence  is  not  as 
significant  in  the  RMS  calculations  from  which  signal  degradation  was  estimated. 

D.  SIGNAL-TO-NOISE  IMPROVEMENT 

1.  Improvement  from  Single  Sensor  to  Subarray  Beam 

In  forming  subarray  beams,  noise  levels  were  reduced  about  7.5  dB, 
slightly  less  than  the  7.  8  dB  expected  from  6  equalized  sensors.  Signals  were  de- 
graded  anyrrhere  from  0  dB  to  2  dB,  with  an  average  degradation  of  about  1  dB.  Thu. 

signal-to-noise  improvement  was  approximately  6.  5  dB,  with  variations  of  about 
1  dB  above  or  below  this  figure. 
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TABLE  IV-2 


SIGNAL  DEGRADATION  FROM  AVERAGE  SUBARRAY  BEAM  TO 
ADJUSTED-DELAY  BEAM 


Event 

Broadband  Degradation 
(dB) 

Degradation  for  Filtered 
Beams  (dB) 

KUR/185/15N 

2.2 

2.9 

WRS/277/10N 

8.6 

8.4 

WRS/295/05N 

7.3 

11.8 

AND/087/08N 

0.9 

2.8 

SIN/207/01N 

3.  0 

3.6 

BLS/210/19N 

— 

4.3 
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2.  Improvement  from  Subarray  Beam  to  Adjusted -Delay  Beam 

For  twenty-five  events,  subarray  signal  powers  over  a  short  gate 
starting  just  before  the  signal  arrival  were  averaged  over  all  operating  subarrays. 
Gate  length  varied  from  6.4  seconds  to  19.2  seconds,  but  was  normally  6.4  seconds. 
Noise  powers  for  the  same  subarrays  over  a  longer  gate  starting  about  75  seconds 
and  ending  about  5  seconds  before  the  signal  arrival  also  were  averaged.  From 
these  figures,  an  average  subarray  beam  signal-to -noise  ratio  was  computed  and 
compared  with  a  similarly  computed  signal-to-noise  ratio  for  the  adjusted-delay 
beam.  Figure  IV-4  is  a  histogram  of  the  signal-to-noise  ratio  improvement  from 
average  subarray  beam  to  adjusted-delay  beam  (where  all  the  beams  have  been 
passed  through  the  standard  detection  filter).  Most  commonly,  signal-to-noise 
improvement  is  about  11  dB.  In  the  case  of  three  events  from  western  Russia  and 
southwestern  Russia,  improvement  is  below  3.  5  dB,  probably  because  adjusted- 
delays  are  inaccurate. 

Signal-to-noise  improvement  of  11  dB  from  subarray  to  adjusted- 
delay  beam  is  somewhat  higher  than  would  be  expected  with  the  3.4  dB  signal 
degradation  figure  derived  from  the  events  KUR/185/15N,  AND/087/08N,  SIN/207/01N 
and  BLS/210/19N.  In  fact,  the  improvement  for  these  events  was  10.  6  dB,  10.4  dB, 

10.  5  dB  and  9.  0  dB,  respectively.  It  appears  that  the  3.  4  dB  signal  degradation 
figure  for  these  events  is  about  1  dB  higher  than  normal. 

Single  sensor  to  subarray  signal-to-noise  ratio  improvement  of 
6.5  dB  together  with  subarray  to  adjusted-delay  beam  improvement  of  11  dB  yields 
a  total  signal-to-noise  ratio  improvement  of  17.  5.  Without  signal  degradation,  the 
single- sensor  to  adjusted-delay  beam  noise  reduction  would  have  given  an  improve¬ 
ment  of  21  dB.  It  appears  that  about  1  dB  was  lost  at  the  subarray  level  and  about 
2.  5  dB  in  forming  the  array  beam  from  the  subarray  beams. 

3.  Comparison  of  Signal-to-Noise  Ratios  for  the  Reference  Subarray 
and  Adjusted-Delay  Beam 

Because  of  the  large  variation  in  signal  amplitudes  from  subarray 
to  subarray,  it  is  useful  to  compare  the  reference  subarray  signal-to-noise  ratio 
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Signal -to-Noise  Ratio  Improvement  (dB) 
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with  that  of  the  adjusted -del ay  beam.  Note  that' the  referenc^e  subarray  is  the 

I  ' 

subarray  with  the  highest  signal-to-noise  ratio.  Figure  rV-5  is  a  histogram,  for 
data  passed  through  the  standard  detection  filter,  of  the  signal-to-noise  ratio  im- 

t 

provement  from  the  reference  subarray  bearA  to  the  adjusted -delay  beam  ffar  26 
events.  In  four  cases,  the  signal-to-noise  ratio  is  higher  at  the  reference  sub-  . 

array.  Two  of  these  events  are  from  western  Russia  .and  two  are  from  southwest- 

;  } 

ern  Russia.  The  remaining  22  events  vary  from  0  to  7.  5  dB  increasp  in  signa,!- 
to-noit  'i  ratio  on  the  adjusted-delay  beam.  Two  peaks  in  improvement,  at  1  dB 
and  5  dB,  are  observed,  of  which  the  second  i,s  the  mokt  comiiion.  If  theisubarray 
with  the  highest  signal-to-noise  ratio  were  used  as  the  array  beam  instead  of'the 
adjusted-delay  beam,  it  is  interesting  that  no  more ‘than  7.  5  dB  in  signal-to-noise 
ratio  would  be  lost.  In  a  few  cases,  the  signal-to-noise  ratio  would  be  increased. 
This  fact,  however,  does  not  mean  that  all  but  one  subarray  could  be  eliminated 
from  the  array  with  signal-to-noise  ratio  loss  no  higher'  than  that  indicated  by  the 
figure,  because  the  subarray  with  the  strongest  signal  changes  with  eVent  location. 

I 

4.  Signal-to-Noise  Ratio  Improvement  bbtained  by  Diversity-Stack  , 

Beamforming  ^ 

1 

Selective  weighting  of  the  subarray  beams  according  to  the  signal- 
to-noise  ratio  on  each  prior  to  the  time -shift -and -sum  process  provides  some 
improvement  in  the  signal-to-noise  ratio  of  the  resultant  b^am.  Such  a  beam,  is 
called  a  diversity-stack  beam.  Figure  IV-6  is  a  histogram  of  the  improvement 
achieved  in  this  way  for  26  events  whose  subarray  jaeams  were  filtered  using  the 

I 

standard  detection  filter.  In  most  cases,  improvement  is  close  to  1'  dB.  The 
two  cases  that  resulted  in  degradation  were  caused  by  inability  to  determine  i 
accurately  the  signal  power  on  the  subarray  beam.  For  both  these  eWnts,  the 
average  subarray  signal-to-noise  ratio  was  less  ‘than  0  dB. 
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SECTION  V 

PRELIMINARY  ESTIMATE  OF  THE  NORSAR  SHORT-PERIOD 
TELESEISMIC  DETECTION  THRESHOLD 

A  preliminary  study  of  the  detectability  of  teleseismic  P-waves 
using  the  NORSAR  short-period  array  has  been  made.  The  adjusted-delay  array 
beam  was  used  for  P-wave  detection.  In  almost  all  cases  where  event  magnitude 
was  below  4.  5,  a  "standard"  detection  filter  was  used  (see  Figure  n-4).  The  fre- 
quency  response  of  this  filter  was  -10  (f  .  2)^  dB,  where  f  is  the  frequency  in 
Hz.  To  obtain  this  filter,  a  filter  was  designed  to  maximize  the  ratio  of  the  low¬ 
est  fUtered  signal  power  to  the  highest  fUtered  noise  power  for  an  ensemble  of 
about  25  signals  and  25  noise  samples.  From  the  filter  derived  in  this  way,  the 
"standard"  detection  filter  was  selected  by  requiring  a  flat  response  near  the  cen¬ 
ter  frequency  and  a  rolloff  rate  approximating  that  of  the  designed  filter.  In  sev¬ 
eral  cases,  detections  were  made  using  the  "standard"  filter  when  a  0.  5  to  5.  0  Hz 
bandpass  filter  faUed  to  work.  There  were  only  one  or  two  cases  where  a  differ¬ 
ent  (narrower  band)  filter  allowed  detection;  thus  the  standard  filter  is  close  to 
optimum  for  detection  of  Eurasian  events. 

The  data  base  used  in  this  study  consisted  of  103  events  reported 
in  the  LASA  and/or  the  PDE  bulletin.  Of  this  population  a  total  of  nineteen  events 
ranging  in  body-wave  magnitude  from  3.  5  to  4.  3  were  not  detected  by  NORSAR. 

The  detected/ not  detected  status  of  a  subset  of  this  population  as  a  function  of  m 
and  distance  from  NORSAR  is  Ulustrated  in  Figure  V-1.  The  expected  tendency'’ 
for  detectability  to  decrease  as  a  function  of  increasing  epicentral  distance  is  not 
clearly  shown  by  this  figure.  This  is  possibly  due  at  least  in  part  to  the  scarcity 
of  data,  but  several  additional  observations  should  be  made.  Four  events  in  the 
epicentral  distance  range  twenty  to  thirty  degrees  and  ranging  in  magnitude  from 
.  8  to  4.  1  were  not  detected.  Three  of  these  events  occurred  in  southwest  Russia. 
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Considerable  difficulty  has  been  encountered  in  forming  beams  for  this  particular 
source  region.  Signal  degradation  by  the  subarray  beams  ranges  as  high  as  four 
dB,  and  additional  degradation  results  from  forming  the  array  beam.  As  a  result 
the  signal-to -noise  ratio  increase  provided  by  array  processing  for  this  source 
region  is  eight  or  nine  dB  lower  than  for  other  source  regions.  Poor  detectability 
then  for  this  region  is  not  surprising. 

The  second  striking  feature  of  Figure  V-1  is  the  mixture  of  detec¬ 
tions  and  non-detections  in  the  distance  range  sixty  to  eighty  degrees.  Based  on 
these  events,  the  detection  probability  does  not  appear  to  vary  as  a  function  of 
magnitude  over  the  range  3.6  to  4.  3.  These  events  are  almost  exclusively  from 
the  Japan-Kuriles -Kamchatka  arc.  The  data  were  examined  in  search  of  a  pattern 
relating  to  precise  epicentral  location,  depth  of  source,  and  background  noise  level. 
No  such  pattern  was  evident.  This  seemingly  anomalous  behavior  is  not  understood 
at  the  present  time. 

In  view  of  these  considerations  it  is  desirable  to  treat  detection 
probability  as  a  function  of  source  region.  Unfortunately,  the  size  of  the  population 
at  the  present  time  does  not  permit  this,  and  it  is  necessary  to  lump  the  entire 
population  together  in  order  to  obtain  a  reasonably  stable  estimate  of  the  detection 
threshold.  As  more  events  are  processed,  it  will  be  possible  to  refine  these  re¬ 
sults  and  obtain  regional  detectability  estimates. 

Figure  V-Z  is  a  histogram  representing  the  magnitude  distribution 
of  the  103  events  which  were  processed.  For  magnitudes  at  which  not  all  processed 
events  were  detected,  the  number  of  events  detected  is  indicated  by  a  dashed 
line.  These  data  were  used  to  obtain  the  estimates  of  incremental  detection  pro¬ 
bability  given  in  Figure  V-3.  In  an  effort  to  stabilize  these  estimates,  the  events 
were  grouped  in  bins  0.  2  magnitude  units  wide.  The  detection  probability  in  each 
bin  is  the  raMo  of  the  number  of  events  detected  to  the  total  number  of  events  pro¬ 
cessed.  As  mentioned  above,  the  population  is  sparse  and  there  is  considerable 
resultant  scatter  in  the  figure.  It  would  appear,  however,  that  the  90%  detection 
level  is  somewhat  near  m^=  4.  3  or  4.  4,  which  is  about  0.  5  m^^  units  above  LASA. 
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INCREMENTAL,  DETECTION  PROBABILITY  OF  THE 


As  discussed  in  Section  II,  about  0.  3  units  can  be  accounted  for  on  the  basis 
of  higher  NORSAR  noise  levels  in  the  detection  band.  The  remaining  0.  2  uiiitd 

i 

may  be  real,  may  be  due  to  the  lack  of  data  or  may  be  associated  with  bias  re¬ 
sulting  from  the  four  close-in  events  discussed  above.  i  i 

i 

! 

It  should  be  noted  that  in  most  cases  the  adjusted  delays  u^ed  in 
beamforming  were  estimated  directly  from  the  event  itself,  xkus  for  large  evehts 
these  adjusted  delays  should  be  the  best  available.  For  the  smaller  events,  how¬ 
ever,  there  exists  the  possibility  of  contamination  of  the  adjusted  de^ys  by  back¬ 
ground  noise.  The  question  as  to  whether  the  use  of  regional  adjusted  delays,  es¬ 
timated  from  an  adequate  population  of  large  events,  will  improve  br  degrade  these 
results  at  NORSAR  has  not  been  fully  answered.  A  second  qualification  of  these 
results  arises  from  the  fact  that  in  all  cases  the  data  were' beamsteered  for  a  known 
source  location  and  the  detection/ no  detection  decision  was  ma^e  by  analyst.  Thus 
the  degradation  inherent  in  an  automated  detection  system  was  not  a  factor  here.  ' 

In  view  of  these  considerations  the  estimate  of  detection  threshold  must  be  consid¬ 
ered  as  only  tentative.  , 

j  1 

Future  work  in  this  area  will  concentrate  on  the  following  gojals: 

•  Expansion  of  the  population  of  events  processed.  ' 

•  Determination  of  regional  detection  thresholds,  and  resolution  of 
anomalies  within  the  regions. 

•  Evaluation  of  the  efficacy  of  regional  adjusted  delays  for  use  in 
NORSAR  beamforming. 

:  i 
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SECTION  VI 

, ,  I  > 

short -PERIOD  DISCRIMINATION 
A.  DEFINITION  OF  DISCRIMINANTS 

'  I 

I  I 

1.  P30  Mean  Square  ' 

'  1 

I  I  ■  : 

This  discriminant;  which  is  a  measure  of  event  complexity,  is  com¬ 
puted  by  crosscorrelat^ng  4  sec  of  the  waveform  (beginning  a  few  points  before  P- 
wave  onset)  with  the  next  30  seconds  ;of  the  waveform  and  with  the  noise  preceding 
the  signal.  A  mean  square,  weighted  by  the  lag,  is  then  computed  from  the  corre- 

i 

'  lations  over  both  30  seconds  pf  the  noise  and  30  seconds  of  the  signal.  The  noise 
mean  square  is  subtracted  from  the  signal  mean  square  to  obtain  the  discrimi¬ 
nant  used  (Texas  Instruments  Incorporated,  1971). 

I  2.  Autocorrelation  Mean  Square 

J 

This  discriminant  is  also  a  measure  of  complexity.  The  autocorrela¬ 
tions  of  a  30-sbcondnoise  gate  and  ofa  30-second  signal  gate  are  computedanda 
weightedmean  square  then  derived  from  these  correlations  for  the  noise  and  signal. 
Th|e  discriminant  is  derived  from  the  signal  mean  square  minus  the  noise  mean. 

'  I 

3.  I  Eiivelope  Difference 

'  , 

This  discriminant  is  also  derived  from  the  P30  correlation  by 

.  I 

computing  the  mean- square  difference  between  the  envelope  correlation  and  a 
fixed  decaying  exponential,  the  decay  rate  of  which  is  the  average  rate  for  an 
ensemble  of  16  explosions  recorded  at  LASA.  As  with  the  first  two  statistics, 
envelope  difference  is  a  measure  of  complexity.  i 

i  '  ' 

* 

4.  Dominant  Period 

I 

i 

This  discriminant  is  computed  by  finding  the  cycle  in  the  waveform 

I  I 

with  the  niaximum  absolute  amplitude;  the  dominant  period  is  the  duration  of  this 

cycle  in  seconds.  This  parameter  can  be  estimated  with  some  confidence,  even 

’1 

j 

!  '  '  ■ 
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for  events  with  a  relatively  low  signal -to -noise  ratio.  The  dominant  period  dis¬ 
criminant  is  a  rough  measure  of  spectral  energy  distribution. 

5.  Spectral  Ratio 

This  discriminant  is  derived  from  the  signal  power  spectrum  over 
a  gate  beginning  just  before  the  signal  arrival.  In  the  original  short-period  dis¬ 
crimination  package,  the  power  spectrum  is  smoothed  over  three  frequency  points, 
and  the  power  in  two  bands  (0.  35  to  0.  85  Hz  and  1.  45  to  1.  85  Hz)  is  computed. 
Dividing  the  power  in  the  first  band  by  that  in  the  other  band  gives  the  discrimi¬ 
nant  used.  The  bands  just  mentioned  appear  to  produce  the  best  separation  of 
earthquakes  and  presumed  explosions  for  events  processed  using  LASA  data. 

Noise  spectra  from  NORSAR  short-period  data,  however,  have 
much  higher  relative  levels  at  frequencies  below  0.5  Hz.  Most  shallow  earth¬ 
quake  spectra,  furthermore,  drop  rapidly  at  frequencies  above  1.5  Hz,  while 
presumed  explosions  do  not  begin  to  fall  off  significantly  until  about  2.  5  Hz 
Thus,  the  NORSAR  spectral  ratio  discriminant  was  computed  by  taking  the  ratio 
of  energy  in  the  band  1.  5  Hz  to  5.  0  Hz  to  that  in  the  band  0.  55  Hz  to  1.5  Hz. 

B.  NORSAR  SHORT -PERIOD  DISCRIMINATION  RESULTS 

Short-period  discriminant  values  fov  the  discriminants  just  de¬ 
fined  are  plotted  as  a  function  of  body-wave  magnitude  for  all  85  detected  events 
in  Figures  VI- 1  to  VI- 12.  Shallow  earthquakes  and  earthquakes  of  unknown  depth 
are  represented  by  a  cross.  Deep  earthquakes  (of  depths  greater  than  100  km) 
are  denoted  by  a  triangle.  Presumed  explosions  are  indicated  by  an  asterisk. 

All  seven  events  from  the  Western  Hemisphere  in  these  graphs  are  surrounded 
by  a  circle.  The  figures  show  that,  as  has  been  previously  observed  (Lacoss  1969), 
the  short  period  discriminants  are  not  effective  for  deep  earthquakes.  The  figures 
also  show  that  the  discriminant  values  are  very  different  for  the  limited  Western 
Hemisphere  earthquake  and  presumed  explosion  ensemble  than  they  are  for  the 
Eurasian  area.  This  is  a  very  important  (and  sometimes  overlooked)  observation; 
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because  the  basic  separation  of  the  two  populations  usually  is  not  large  (for  a 
given  SP  discriminant)  path  effects  can  "overpower"  any  separation  inherent 
in  the  discriminant.  Thus,  a  SP  discriminant  which  works  well  at  a  station  for 
source  area  A  may  not  work  at  all  for  source  area  B.  Conversely,  an  SP  dis¬ 
criminant  which  works  well  for  a  source  area  at  station  X  may  not  work  at  all 
for  station  Y. 

The  remainder  of  this  section  is  confined  to  a  discussion  of  the 
performance  of  the  SP  discriminants  for  shallow  Eurasian  events  (crosses  in 
the  figures)  and  presumed  explosions  (asterisks  in  the  figures). 

Figure  VI  1  shows  the  logarithm  of  the  P30  mean  square  discrimi¬ 
nant  as  a  function  of  body-wave  magnitude  for  the  reference  subarray  beam;  there 
is  a  wide  region  of  overlap  between  the  earthquake  and  presumed  explosion 
population.  Figure  VI-2  shows  the  same  discriminant  for  the  adjusted”delay  beam 
separation  is  somewhat  better,  but  thei'e  is  still  a  considerable  region  of  overlap. 

Figure  VI-3  shows  the  log  of  the  autocorrelation  mean  square  dis¬ 
criminant  as  a  function  of  body-wave  magnitude  for  the  reference  subarray  beam. 
Again,  there  is  a  wide  region  of  overlap  between  the  earthquakes  and  presumed 
explosions.  When  the  same  discriminant  is  plotted  for  the  adjusted-delay  beam 
(Figure  VI-4),  better  separation  is  apparent,  but  overlap  still  occurs. 

The  last  of  the  discriminants  based  on  event  complexity,  envelope 
difference,  is  presented  in  Figure  VI-5  for  the  reference  subarray  beam.  Over¬ 
lap  still  persists.  In  Figure  VI-6  the  same  discriminant  computed  from  the 
adjusted -delay  beam  is  given.  Slightly  better  separation  is  apparent,  but  overlap 
still  is  severe. 

The  remaining  discriminants  to  be  presented  in  this  section  are 
based  on  spectral  energy  distribution.  The  first  of  these,  the  dominant-period 
discriminant,  is  shown  in  Figure  VI-7.  It  would  be  possible  to  draw  a  straight 
line  separating  all  but  one  earthquake  from  the  presumed  explosions.  But  the 
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Magnitude 
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FIGURE  VI- 1 
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FIGURE  VI -2 

ADJUSTED-DELAY  BEAM  P30  MEAN  SQUARE  DISCRIMINANT 
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FIGURE  VI-3 
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FIGURE  VI-5 
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FIGURE  VI-6 
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FIGURE  VI-7 

REFERENCE  DOMINANT  PERIOD  (SECONDS) 


separation,  even  ignoring  the  one  earthquake,  is  extremely  small.  If  presumed 
explosions  of  smaller  magnitude  had  been  processed,  it  is  likely  that  there  would 
be  more  overlap.  The  dominant-period  discriminant  calculated  from  the  adjusted- 
delay  beam  is  presented  in  Figure  VI-8.  In  this  case,  there  is  one  presumed 
explosion  inside  the  earthquake  population.  The  reason  for  this  anomalous  point 
will  be  discussed  later.  Ignoring  this  point,  there  is  a  noticeable  separation 
between  the  earthquakes  and  the  presumed  explosions. 

The  spectral  ratio  of  energy  in  the  band  0.  3  to  0.  8  Hz  to  that  in 
the  band  1.  4  to  1.8  Hz  is  graphed  in  Figure  VI-9  for  the  reference  subarray  beam. 
The  choice  of  bands  was  based  on  events  processed  using  LASA  data.  Two  pre¬ 
sumed  explosions  lie  within  the  earthquake  population,  and  one  earthquake  lies 
within  the  presumed  explosion  population.  In  Figure  VI- 10,  the  same  discrimi¬ 
nant  is  plotted  for  the  adjusted-delay  beam.  Here  the  overlap  is  considerably 
worse. 

After  examination  of  spectra  for  earthquakes  and  presumed 
explosions  received  at  NORSAR,  the  bands  1.  5  to  5.  0  Hz  and  0.  55  to  1.  5  Hz  were 
selected  for  use  in  a  spectral  ratio  discriminant.  Figure  VI- 11  is  a  plot  of  the 
resultant  spectral  ratios  for  the  reference  subarray  beam.  Two  presumed  explo¬ 
sions  lie  within  the  earthquake  population.  The  corresponding  spectral  ratios 
for  the  adjusted-delay  beam  are  shown  in  Figure  VI- 12.  The  same  two  presumed 
explosions  lie  within  the  earthquake  population.  In  both  figures,  all  other  pre¬ 
sumed  explosions  from  the  Eurasian  area  lie  along  a  curve  at  the  top  right  of  the 
graphs. 

One  of  the  two  presumed  explosions  within  the  earthquake  popula¬ 
tion  was  a  magnitude  5.3  event  from  Western  Russia  on  day  295  of  1971 
(Figure  VI- 13).  It  was  very  difficult  to  decide  on  the  adjusted  delays  to  be  used 
in  forming  the  array  beam  for  this  event.  When  different  portions  of  the  signal 
trace  or  different  features  in  the  trace  were  used  for  determining  the  adjusted 
delays,  different  results  were  obtained.  It  is  likely  that  much  of  the  high-frequency 
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FIGURE  VI- 10 
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adjusted-delay  beam  spectral  ratio 
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energy  was  lost  in  forming  the  array  beam.  In  Figure  VI-7  (dominant  period  for 
the  reference  subarray),  this  event  lies  at  the  edge  of  the  presumed  explosion 
population.  In  Figure  VI-8  (dominant  period  for  the  adjusted-delay  beam),  how- 

i 

ever,  this  event  is  the  one  anomalous  point.  v 

The  other  presumed  explosion  within  the  earthquake  population 
was  a  magnitude  5.4  event  from  Eastern  Kazakh  on. day  181  of  1971  (see  Figure 
VI- 14).  The  dominant-period  discriminant  is  just  able  to  place  this;event  within 
the  presumed  explosion  population  (see  Figures  VI-7  and  VI-8),  whereas  the 
spectral  ratio  discriminant  places  this  event  within  the  earthquake  population. 

A  comparison  of  the  reference  subarray  beam  and  the  adjusted-delay  beam 
(Figure  VI- 14)  would  suggest  that  the  beamforming  process  was  responsible  for 
the  loss  of  high-frequency  energy.  In  the  case  of  the  reference  subarray  beam, 
the  spectral  ratio  for  this  event  (see  Figure  VI-11)  lies  just  within  the  earth-  * 
quake  population.  For  the  adjusted-delay  beam,  however,  the  spectral  ratio 
for  this  event  (see  Figure  VI- 12)  is  buried  somewhat  nore  deeply  within  the 
earthquake  population.  The  adjusted  delays  for  this  area  are  believed  to  be 
well -determined.  We  cannot  explain  the  reason  why  so  much  of  the  high- 
frequency  energy  is  lost  in  the  adjusted-delay  beam. 
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SECTION  vn 

CONCLUSIONS  AND  FUTURE  PLANS 


CONCLUSIONS 


Conclusions  about  the  performance  of  the  short  period  NORSAR 

array,  based  on  analysis  of  just  over  100  signals  (primarily  from  Eurasia)  and 
72  noise  samples,  are  given  below. 


Data  quality  is  excellent.  For  about  one-haU  of  the  samples  all 
sensors  were  operational.  For  most  of  the  rest  of  the  samples  one  or  two 
aubarrays  (6  or  12  sensors)  were  inoperative;  the  worst  data  loss  was  19  sensors. 
There  were  essentially  no  spikes  or  glitches  in  the  data. 

Conclusions  from  the  noise  analysis  are  based  on  72  samples  and 
can  be  stated  with  good  confidence.  Major  results  are: 

The  noise  spectral  shape  is  very  simple.  The  peak  occurs  at  3  to 

6  seconds  and  the  spectra  decrease  rapidly  at  shorter  periods.  The 

spectral  shape  does  not  change  significantly  either  across  the  array 
or  with  time. 

.  Noise  levels  are  very  similar  across  the  array.  Maximum  single 

sensor  variations  typically  are  ±  6  dB,  and  most  sensors  are  with- 
in  ±  3  dB  of  the  average  single  sensor  level.  Variation  am-,ng  sub¬ 
array  beam  noise  levels  in  the  region  of  significant  power  (0  to  2.  0 
Hz)  is  i  2  dB. 

Wideband  RMS  noise  levels  show  a  definite  seasonal  dependence- 
wintertime  levels  are  6  dB  higher.  The  increased  wintertime 
level  results  because  the  3  to  6  second  microseismic  peak  is 
stronger  in  the  winter.  Application  of  the  "standard"  bandpass  filter, 
which  roUs  off  sharply  at  high  frequencies,  reduces  the  difference 
between  summertime  and  wintertime  noise  levels. 
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•  RMS  levels  on  the  "standard">filtered  adjusted-delay  array  beam 
range  from  0.  04  to  0.  17  m^,  and  typically  are  0.  12  m^.  This  is 
about  a  factor  of  two  higher  than  LASA  detection-filtered  beam 
noise  levels. 

•  Multiple  coherence  levels  witbin  a  subarray  are  low  except  at  the 
3  to  6  second  microseismic  peak.  Inter -subarray  multiple  coher¬ 
encies  are  low  over  the  entire  0  to  5.  0  Hz  band. 

Major  conclusions  from  the  signal  analysis  are: 

•  Except  for  a  few  close-in,  high-frequency  events  signal  similarity 
is  good  within  a  subarray.  Among  subarrays,  however,  similarity 
is  quite  variable. 

•  Amplitude  variations  across  the  array  typically  are  4;1.  Some 
regions  (e.  g.  ,  Kazakh)  show  variations  as  high  as  10:1. 

•  Eurasian  signals  usuedly  have  a  substantial  amount  of  high  fre¬ 
quency  energy  (out  to  about  2  Hz).  Spectral  shapes  are  quite  vari¬ 
able, but  the  close-in  events  (  A <30°)  generadly  have  more  high- 
frequency  content.  The  limited  ensemble  of  Western  hemisphere 
events  processed  show  substantially  less  high  frequency  energy 
than  the  Eurasian  events. 

•  Time-delay  anomalies  (deviation  from  plane  wave  propagation  along 
the  great  circle  path)  are  not  significant  for  subarray  beamforming, 
except  perhaps  for  a  few  close-in  events.  Anomalies  are  significant 
however,  among  subarrays  and  corrections  are  essential  for  array 
beamforming.  Consistent  subarray  delay  anomalies  could  be  ob¬ 
tained  for  all  regions  except  those  within  30°  epicentral  distance  of 
NORSAR. 

•  NORSAR  m^'s  averaged  about  0.2  units  less  than  PDE  (and  LASA) 
m^*s.  This  discrepancy  can  be  explained  as  signal  loss  in  array 
beamforming,  and  thus  it  appears  that  NORSAR  m^'s,  if  corrected 
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for  signal  loss,  are  about  the  same  as  PDE  The  PDE- 

NORSAR  m^  differences  do  appear  to  be  larger  at  low  magnitudes; 
this  may  be  because  the  low  magnitude  PDE  ni^'s  are  obtained  from 
the  quieter  stations  in  the  PDE  network  (e.  g. ,  TFO),  which  normally 
have  lower  thcin  average  signal  amplitudes. 

Major  conclusions  from  the  array  processing  performance  study 

are: 

•  VFT  noise  rejection  is  achieved  over  the  entire  0  to  5  Hz  band  for 
both  subarray  and  array  beamforming.  Thus  noise  rejection  totals 
about  21  dB  (8  for  subarray  and  13  for  array  beamforming). 

•  Signal  degradation  for  subarray  beamforming  typically  is  1  dB  but 
for  some  close-in  high-frequency  signals  loss  appears  to  be  about 
3  or  4  dB. 

•  Signal  degradation  for  arr^y  beamforming  is  quite  variable,  but 
generally  appears  to  be  only  about  2  to  4  dB.  On  some  events, 
however,  degradation  appears  to  be  as  high  as  11  dB.  It  should 

be  pointed  out  that  it  is  quite  difficult  to  estimate  array  beamform¬ 
ing  signal  degradation,eBp<;cially  when  it  is  large,  because  of  the 
dissimilarity  of  subarray  waveforms. 

•  Net  signal-to-noise  ratio  (SNR)  typically  is  15  to  18  dB.  For 

the  problem  close-in  events  gains  drop  as  low  as  about  10  dB.  For 
these  events  the  reference  subarray  beam  often  has  a  higher  SNR 
than  the  array  beam.  (Usually  the  array  beam  SNR  is  about  1  to  5 
dB  higher  than  the  reference  subarray  beam  SNR). 

•  Diversity  stack  beamforming  provides  0  to  2.  5  dB,  and  typically 
1  dB,  better  SNR  improvement  than  the  adjusted -delay  beam. 

•  For  detection  of  Eurasian  events,  a  bandpass  filter  with  corner 
frequencies  at  about  1.2  and  2.8  Hz  and  a  very  sharp  rolloff  at  low 
frequencies  appears  to  be  about  optimum.  The  relatively  high  band¬ 
pass  is  desirable  because  event  SNR's  generally  peak  at  about  1.5  Hz. 
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This  filter  may  not  be  best  for  Western  hemisphere  events  which 
appear  to  have  somewhat  lower  signal  spectral  content. 

A  preliminary  estimate  of  the  NORSAR  detection  threshold  for 
Eurasian  events  gives  a  90%  incremental  value  of  m^-  4.  3.  Much  more  data 
are  required  to  obtain  a  reliable  estimate,  however. 

Standard  short-period  discriminants  were  calculated  for  all  detec¬ 
ted  events,  and  the  conclusions  are: 

•  Discriminants  based  on  event  complexity  (P30  mean  square,  auto¬ 
correlation  mean  square,  and  envelope  difference)  do  not  appear 
to  be  very  effective  in  separating  shallow  earthquakes  or  earth¬ 
quakes  of  unknown  depth  from  presumed  explosions  for  Eurasian 
events. 

•  Discriminants  based  on  spectral  energy  distribution  (dominant 
period  and  spectral  ratio)  appear  to  work  reasonably  effectively, 
although  complete  separation  between  Eurasian  earthquakes  and 
presumed  explosions  was  not  achieved.  Surprisingly,  the  very 
simple  dominant  period  measurement  gave  the  best  results.  The 
spectral  ratio  of  energy  in  the  bands  0.  3  to  0,  8  Hz  and  1.4  to  1.8 
Hz  was  clearly  inferior  to  the  ratio  of  energy  in  the  bands  0.  55  to 
1,5  Hz  and  1,5  to  5,  0  Hz,  This  fact  is  not  surprising,  since  the 
first  two  bands  were  based  on  events  processed  using  LASA  data. 

•  The  short  period  discriminant  values  for  Eurasian  events  are  sign¬ 
ificantly  different  that  those  obtained  for  a  limited  ensemble  of 
Western  hemisphere  earthquakes  and  presumed  explosions.  This 
observation  points  out  that  the  effectiveness  of  short-period  dis¬ 
criminants  depends  on  the  source  and  station  location. 
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B. 


FUTURE  PLANS 


Future  NORSAR  short-period  evaluation  efforts  will  concentrate 
on  increasing  the  Eurasian  event  ensemble,  emphasizing  low-magnitude  events, 
in  order  to  improve  estimates  of  the  array  detection  capability.  Close-in  events 
will  be  analyzed  in  detail  in  an  attempt  to  find  techniques  for  improving  array 
SNR  gain. 
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